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Abstract. We report here on recent progress in understanding the birth conditions of neutron stars and the way how 
supernovae explode. More sophisticated numerical models have led to the discovery of new phenomena in the supernova core, 
for example a generic hydrodynamic instability of the stagnant supernova shock against low-mode nonradial deformation and 
the excitation of gravity- wave activity in the surface and core of the nascent neutron star. Both can have supportive or decisive 
influence on the inauguration of the explosion, the former by improving the conditions for energy deposition by neutrino 
heating in the postshock gas, the latter by supplying the developing blast with a flux of acoustic power that adds to the energy 
transfer by neutrinos. While recent two-dimensional models suggest that the neutrino-driven mechanism may be viable for 
stars from ^SM© to at least ISM©, acoustic energy input has been advocated as an alternative if neutrino heating fails. 
Magnetohydrodynamic effects constitute another way to trigger explosions in connection with the collapse of sufficiently 
rapidly rotating stellar cores, perhaps linked to the birth of magnetars. The global explosion asymmetries seen in the recent 
simulations offer an explanation of even the highest measured kick velocities of young neutron stars. 
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INTRODUCTION 

Improved numerical tools and the increasing power 
of modern supercomputers have brought considerable 
progress in modeling stellar core collapse in the past 
years. It is possible now to simulate the complex phys- 
ical processes in the deep interior of supernovae with un- 
precedented sophistication and detailedness. 

It has become clear meanwhile that the explosions of 
massive stars are a generically multi-dimensional phe- 
nomenon. This insight was fostered by the fact that 
spherically symmetric (ID) simulations, which became 
available with a fully energy dependent solution of the 
Boltzmann transport problem for the neutrinos only re- 
cently (see 1 1] for an overview and comparison of differ- 
ent numerical approaches) confirmed and solidified older 
ID results of the 1980's and 1990's, namely that explo- 
sions in the ID models could not be obtained, neither by 
the prompt bounce- shock nor by the delayed neutrino- 
heating mechanism, at least not for progenitor stars of 
more than lOM© and on a timescale of roughly one sec- 
ond after core bounce 10, S 0, [sl]. Moreover, the lat- 
est generation of multi-dimensional simulations has pro- 
vided evidence for a variety of routes that can lead to ex- 
plosions when nonradial phenomena are accounted for. 
These routes seem to depend on the properties and condi- 
tions in the progenitor stars like their mass and structure 
and the amount of angular momentum in their core. 

In the following we will briefly review these recent de- 
velopments in the multi-dimensional modeling of stellar 
core collapse and explosion, and we will critically dis- 
cuss the status of the present simulations. 



A better understanding of the explosion mechanism of 
core-collapse supernovae is not only important for inter- 
preting the observable properties of the blast, for pre- 
dicting gravitational- wave and neutrino signals, and for 
determining the conditions of nucleosynthesis processes 
that occur during the explosion. It is also and in particular 
essential for establishing the link between the progenitor 
stars and their compact remnants, thus answering ques- 
tions like that of the mass distribution of neutron stars 
and of the stellar mass limit for black hole formation, 
which may happen either directly during the core col- 
lapse or by later massive fallback when the disrupted star 
does not become completely unbound during the explo- 
sion. So far, estimates for such scenarios have been made 
only on the basis of still rather crude self-consistent ex- 
plosion simulations | 6] or by invoking assumptions about 
the mass cut and energy in models with piston-driven ar- 
tificial explosions (e.g., f^). 

BRIEF HISTORICAL EXCURSION 

Due to the huge gravitational binding energy liberated 
in neutrinos, which carry away hundred times more en- 
ergy than needed for the explosion, these particles have 
long been speculated to be the driving agent of the stel- 
lar explosion. Colgate and White 1 8] in a seminal paper 
in 1966 not only proposed gravitational binding energy 
to be the primary energy source of core-collapse super- 
novae, but also that the intense flux of escaping neutri- 
nos transfers the energy from the imploding core to the 
ejected stellar mantle. Nearly twenty years later, Bethe 



and Wilson ^ were the first who described in detail 
the way how this might happen, interpeting thereby the 
physics that played a role in hydrodynamic simulations 
performed by Wilson. They concluded that electron neu- 
trino and antineutrino absorptions on the free neutrons 
and protons that are abundantly present in the shock- 
dissociated matter behind the stalled accretion shock in 
the supernova core, are the primary agents of the energy 
transfer. 

These pioneering computer simulations of the so- 
called delayed neutrino-driven explosion mechanism 
were still conducted in spherical symmetry. The mech- 
anism turned out to be successful only when the neu- 
tron star was assumed to become a more luminous neu- 
trino source by mixing instabilities accelerating the en- 
ergy transport out of its dense interior. The thus enhanced 
neutrino emission led to stronger neutrino heating in the 
overlying layers of the exploding star. Theoretical studies 
and multi-dimensional computer models, however, sug- 
gest that convection and mixing instabilities inside the 
neutron star do not have the necessary big effect (see, e.g. 

Instead, the first multi-dimensional sim- 
ulations, which became available only in the mid 1990's, 
demonstrated that the neutrino-heated layer around the 
forming neutron star is unstable to vigorous convective 
overturn iH HI [H, B E [TtI] • This can raise the ef- 
ficiency of the neutrino energy deposition and thus can 
have a supportive effect on the supernova explosion. The 
first such two-dimensional (i.e. axisymmetric) and three- 
dimensional simulations, however, suffered from a se- 
vere drawback: the physics of the neutrino transport and 
of neutrino-matter interactions, which is essential for dis- 
cussing the power input to the explosion, is so complex 
that it could be treated only in a grossly simplified way. 
In the best models at that time this was done by the so- 
called "grey diffusion approximation". This means that 
the energy-dependence of the neutrino interactions (the 
cross sections of the most important neutrino processes 
typically scale with the squared neutrino energy) was 
ignored and replaced by a "grey" (spectrally averaged) 
description. Moreover, the spatial propagation was ap- 
proximated by assuming that neutrinos diffuse through 
the dense neutron star medium until they decouple and 
stream away from a chosen position, usually from a layer 
somewhat outside of the "neutrino sphere", close to the 
surface of the compact remnant. The historical develop- 
ment of these theoretical studies of the supernova explo- 
sion mechanism is resumed in a recent review 1 18]. 

RECENT RESULTS 

Only in the past years the neutrino treatment in multi- 
dimensional hydrodynamic and magnetohydrodynamic 
(MHD) models of supernovae has seen significant im- 



provements. However, a rigorous solution of the Boltz- 
mann transport equation is still much too time consum- 
ing to be applied in full-scale simulations. Even in axi- 
ally symmetric (2D) models the transport poses a five- 
dimensional problem (see, e.g., l20]), in three dimen- 
sional hydrodynamic simulations it would constitute a 
time-dependent six-dimensional problem. 

All active groups therefore still have to accept some 
simplifications and nevertheless the neutrino transport 
module dominates the computing time for supernova 
simulations by far. The approximations taken by dif- 
ferent groups differ significantly. While the Tucson- 
Jerusalem collaboration employs a 2D flux-limited dif- 
fusion scheme and treats the neutrino energy groups un- 
coupled (an approach that is known from ID simulations 
to be unable to capture important physics), thus gaining 
a modest amount of straightforward parallelism for their 
computations (e.g., |21, 22, 23], the transport treatment 
of the Garching group accounts for the full energy depen- 
dence of the problem and solves on each angular (lateral) 
bin of the 2D grid a full one-dimensional transport prob- 
lem by iterating the moment equations of neutrino num- 
ber, energy, and momentum with a variable Eddington 
factor for the closure that is obtained from the solution of 
a model Boltzmann equation. Moreover, neutrino pres- 
sure gradients and advection in the lateral direction are 
included in this so-called "ray-by-ray plus" approxima- 
tion f2?, . While this approach assumes that the neu- 
trino flux components in lateral direction are zero (i.e., 
the neutrino intensity is taken to be symmetric around 
the radial direction), it allows one to properly treat the 
gradual transition of neutrinos from diffusion at the high 
densities in the neutron star interior to free streaming in 
the much more dilute stellar layers far outside of the neu- 
tron star. This approach has also the big advantage of be- 
ing a direct generalization of the ID case and therefore 
it enables a detailed, well constrained comparison of ID 
and 2D simulations. 

In the following we will summarize the essentials of 
recent two-dimensional studies that have made use of the 
mentioned improvements in the neutrino transport, and 
which have thus contributed to a better understanding of 
the question how the collapse of stellar cores could be re- 
versed to an explosion. The basic requirement for this to 
happen is that some energy reservoir that takes up grav- 
itational binding energy released during stellar core col- 
lapse can be effectively tapped and transferred to matter 
that can get expelled in the explosion. This can happen 
by neutrinos in the context of the neutrino-heating mech- 
anism, but it can also be achieved by magnetic fields in 
magnetohydrodynamic (MHD) explosions. Or it may oc- 
cur, as recently proposed llilll^], through sound waves 
created by violently turbulent gas motions around the im- 
pact sites of accretion downflows on the neutron star sur- 
face, or even by large- amplitude g-mode pulsations of 
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FIGURE 1. Snaphots showing the gas entropy (left half of panels) and the electron-to-nucleon ratio (right half of panels) for the 
explosion of an ^9Mq star with an 0-Ne-Mg core. The plots correspond (from top left to bottom right) to times of 0.097, 0. 144, and 
0.262 seconds after the launch of the supernova shock front and the onset of neutron star formation at the moment of core bounce. 
Note the different radial and color scales of the four panels. Due to the rapid expansion of the shock and of the shock- accelerated 
ejecta into the extremely dilute layers surrounding the 0-Ne-Mg core, the convective pattern freezes out quickly and begins a nearly 
self-similar expansion. The characteristic wavelength of convective structures is roughly 30-45 degrees (corresponding to dominant 
spherical harmonics modes of / = 4,5) and there is no strong contribution of dipolar and quadrupolar asymmetries. 




FIGURE 2. Radii of the supernova shock as functions of 
time for one- and two-dimensional simulations (red and black 
lines, respectively) of the explosion of a star with 0-Ne-Mg 
core. Note that the progenitor used in the 2D simulation was an 
8.8M0 model with an artificially constructed low-density He- 
shell at p < 10^ gcm~^ 1 19], while the ID simulation was per- 
formed with a recently updated progenitor structure in which a 
H-envelope with a much lower density and steeper density de- 
cline was added around the 0-Ne-Mg core (K. Nomoto, private 
communication). This explains the stronger acceleration of the 
shock in the region outside of about 1100 km. 

the neutron star core, leading to acoustically powered ex- 
plosions. 

Neutrino-driven explosions 

Neutrinos extract energy from the huge reservoir of 
degeneracy and thermal energy that is built up inside 
of the nascent neutron star during stellar core collapse. 
These neutrinos diffuse out of the dense interior and be- 
fore streaming off to low-density regions, mostly neutri- 



nos of the electron flavor deposit roughly 10% of their 
energy in the so-called gain layer between the gain ra- 
dius and the stalled supernova shock. Detailed and accu- 
rate numerical models are indispensable to determine the 
exact efficiency of this energy transfer. 

Results obtained by the Garching group for progeni- 
tor stars between 9 and 15 solar masses confirm the vi- 
ability of the neutrino-heating mechanism for triggering 
supernova explosions. However, the inauguration of the 
explosion happens in a different way than expected from 
previous models and the blast properties turn out to dif- 
fer significantly from older calculations with more sim- 
plified neutrino physics. 

Supernova progenitors with less than about lOM© 
clearly vary in their structure and explosion behavior 
from more massive stars. The former class of stars de- 
velops a core composed of oxygen, neon, and magne- 
sium, not of iron, with an extremely steep density gradi- 
ent at its surface. This allows the supernova shock front, 
which is launched at the moment when the neutron star 
begins to form at the center of the collapsing stellar core, 
to expand continuously as it propagates into rapidly di- 
luting inf ailing material. Behind the shock the velocities 
are initially negative, so that no prompt explosion oc- 
curs. But then neutrino heating deposits the energy that 
powers the ensuing blast. Convective overturn develops 
in the neutrino-heated layer behind the outgoing shock 
and imprints inhomogeneities on the ejecta, in entropy 
as well as composition (Fig.[T]). But because of the rapid 
acceleration of the supernova shock and of the postshock 
layer (Fig. O, the pattern of Rayleigh-Taylor structures 
freezes out quickly (when the expansion timescale be- 
comes shorter than the overturn timescale) and the inho- 
mogeneous shell behind the shock begins an essentially 
self- similar expansion. The short convective phase leads 
to large-scale explosion asymmetries, however, without 
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FIGURE 3. Snapshots of the gas entropy for the explosion of a star with 11.2M0 at times 0.14, 0.20, and 0.28 seconds after the 
launch of the supernova shock at core bounce (top left to bottom right). The explosion develops a large bipolar asymmetry although 
the star is not rotating. Note the different radial and entropy scales of the four panels. The first two plots show results from Ref. yfl, 
the last one is from a recent continuation of the same simulation to later times. 
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FIGURE 4. ]V[aximum, average, and minimum radial posi- 
tion of the supernova shock front as functions of time for the 
explosion of the 11.2M0 model displayed in Fig. [3] Note the 
clear signature of several large- amplitude bipolar shock oscil- 
lations due to the standing accretion shock instability (SASI) 
before the blast takes off with an extreme 3:1 deformation. 

global dipolar or quadrupolar deformation (Fig.[T]). 

Evolved stars above roughly 10 Mq produce iron cores 
with a much more shallow density decline outside. Run- 
ning into this denser material damps the initial expan- 
sion of the shock. Moreover, severe energy losses and 
the high mass infall rates cause the shock to even stall at 
a relatively small radius of only about 100 km. Because 
of the small shock stagnation radius, the infall veloci- 
ties of the collapsing stellar core ahead and behind the 
shock are very large. Different from previous calcula- 
tions with simple grey neutrino diffusion, our more so- 
phisticated models show that convection is strongly sup- 
pressed in the rapidly infalling matter behind the shock. 
Neutrino heating is not powerful enough to allow high- 
entropy matter to become buoyant against the accretion 
flow (see 0,0]). Convective overturn behind the shock 
therefore cannot become sufficiently strong to help push- 



ing the stagnant shock farther out and thus to establish 
more favorable conditions for neutrino heating. 

Instead, another kind of nonradial hydrodynamic in- 
stability, the so-called standing accretion shock instabil- 
ity ("SASF'; 1 28]), which can grow efficiently even when 
convection stays weak [29, 30, 27], obtains decisive in- 
fluence on the shock evolution. With highest growth rates 
of the dipole and quadrupole modes |31], it leads to vi- 
olent bipolar sloshing motions of the shock. This drives 
the shock front to larger radii and thus reduces the ac- 
cretion velocities in the postshock layer. The oblique- 
ness of the shock surface relative to the infalling stellar 
core material deflects the accretion flow and stretches its 
path through the layer of neutrino heating. Moreover, the 
SASI also causes secondary convection due to steep en- 
tropy gradients produced in the postshock layer by the 
quasi-periodic expansion and contraction phases of the 
shock. The influence of the SASI thus improves the con- 
ditions for efficient energy deposition by neutrinos, be- 
cause the gas accreted through the stalled shock can stay 
longer in the heating layer and is therefore able to absorb 
more energy from the intense neutrino flux radiated by 
the nascent neutron star 1I27I1 . 

The presence of strong SASI oscillations is visible in 
simulations for an II.2M0 star (for details, see \3]) and 
for a 15M0 star (details in Ref. (H). The SASI turns 
out to be crucial for the explosion in both cases (Figs. [3]- 
[5]). Different from the ^9Mq star, where convection im- 
poses a high-mode asymmetry pattern on the ejecta (see 
Fig.[T]), the preferred growth of the dipole and quadrupole 
(/ = 1 , 2; m = in terms of an expansion of the character- 
istic flow properties in spherical harmonics) modes of the 
SASI leads to a large global anisotropy of the beginning 
supernova blast in the more massive progenitors. 

Our simulations have thus demonstrated the decisive 
role of the standing accretion shock instability in com- 
bination with neutrino heating for initiating neutrino- 
powered supernovae. The explosion may set in with a 
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FIGURE 5. Left: Shock positions near the north pole and near the south pole as functions of time for a 15 Mq star that evolves 
towards the onset of an explosion at > 600ms after core bounce. The gas entropy is color coded. The plot shows many cycles of 
quasi-periodic bipolar shock oscillations due to the standing accretion shock instability (SASI). Right: Gas entropy distribution 
for the \5Mq star at the onset of the explosion at a time of 0.61 s after bounce. A large north-south asymmetry signals a strong 
contribution from the dipole mode. The white line marks the gain radius as lower boundary of the neutrino-heating region (both 
plots are from 1321 ). 



TABLE 1. Estimated explosion properties for 
different progenitor stars with ZAMS mass Mprog 



Mprog 


[ms] 


[Mq] 


[B] 


[Mq] 


^9 


120 


0.01-0.02 


0.2-0.3 


1.36 


11.2 


220 


0.02-0.04 


0.3-0.6 


1.30 


15 


620 


0.05-0.06 


-1.0 


1.55 



* Time of onset of explosion, determined as the moment 
when the total energy in the gain layer, integrated over 
the mass elements with positive specific total energy, 
exceeds 10^^ erg 

Mass of accreted, neutrino-heated, and then ejected 
matter accounting for the explosion energy estimate of 
Eq.O 

** Explosion energy estimated according to Eq. ([T) in 
IB = Ibethe = 10^^ erg 

* Baryonic mass of the neutron star at the onset of the 
explosion 



significant delay after the neutron star begins to form 
(similar results, however with Newtonian gravity instead 
of a relativistic gravitational potential were recently re- 
ported in Ref. 11331] ). For the most massive of the three 
investigated stars this happens about 0.6 seconds later. 
This is not only much later than expected from previous 
simulations, but also constitutes a major computational 
challenge for our 2D modeling with sophisticated multi- 
energy-group neutrino transport, for which reason our set 
of computed cases is still constrained to only three pro- 
genitors. 

Table [T] summarizes the prediced explosion and rem- 
nant properties for these three progenitors. Only in the 
case of the rapidly developing blast of the ~9M0 star 



could the ID and 2D calculations be carried on for a suf- 
ficiently long time to see the explosion energy asymptote. 
For the other two models the explosion energy is built up 
during a possibly long-lasting phase of anisotropic ac- 
cretion and simultaneous expansion of neutrino-heated 
matter, which is a manifestation of the multi-dimensional 
nature of the explosion Q 111 Q HI] . The explosion 
energy can then be roughly estimated from (for details, 
see 1^2]) 

^exp ~ Ey Tacc ~ ^acc , (1) 

where Ey is the net (i.e., heating minus cooling) neutrino 
energy transfer rate to accreted and ejected gas, 

^v-C^acc^v -2X10^1^, (2) 

s 

when Mace ~ O.2M0/S is the mass accretion rate through 
the shock when the explosion begins in the 11.2 and 
15M0 models, ~ 0.5 is the fraction of the accreted 
gas that gets ejected again after being heated by neutri- 
nos, ey lOlMeV/nucleon ~ lO^^erg/g is the average 
specific energy deposited by neutrinos in the gas, and 
Mace = C^acc Tacc is the accreted and subsequently ex- 
pelled gas mass. Accretion can continue until the post- 
shock matter is accelerated to escape velocity by the out- 
going shock, which leads to an estimate of the accretion 
time as Tacc ^ 0-5sMns,i.5V^^9, when Mns,i.5 is the neu- 
tron star mass normalized to 1 .5Mq and Vsh,9 is the shock 
velocity in units of 10^ cm/s (see [32]). 

The non-monotonic behavior of the initial neutron star 
(baryonic) mass is linked to the core size of the progen- 
itors and the surrounding density structure, because the 
latter determines the delay of the explosion and the dura- 
tion of the accretion phase after bounce. 




FIGURE 6. Upper panels: Two snapshots of the entropy distribution at 0.25 s and 1.0 s after bounce for one of the simulations in 
Ref. 1 34]. The anisotropic ejecta distribution leads to a neutron star kick as explained in the text. Lower panels: Neutron star recoil 
velocity (left) and acceleration (right) as functions of time after bounce for the simulation shown in the upper two panels. The solid 
curves are the result deduced directly from the hydrodynamic simulation, the dashed curve in the left plot and the long-dashed curve 
in the right plot are from an independent post-processing analysis of all accelerating effects on the neutron star. It turns out that the 
gravitational pull of the anisotropic ejecta is the main mediator of the neutron star acceleration at ? > 0.5 s (dotted curve in the right 
panel) and dominates anisotropic accretion and outflows and pressure forces (short-dashed, dash-dotted, and dashed-triple dotted, 
respectively, in the right panel) by far. Anisotropic neutrino emission has usually a small influence (difference between solid and 
dotted lines in the left panel). 



Explosion asymmetries and pulsar kicks 

One of the consequences of this SASI- supported 
neutrino-driven mechanism is a global asymmetry of the 
accelerating shock front and of the ejected gas even in the 
absence of rotation (or with only very little rotation) in 
the stellar core. When the dipole and quadrupole modes 
are very strong, the onset of the explosion can resemble 
even a bipolar jet-like blast with a sizable pole-to-equator 
deformation (see Figs. [3] and |4]). 

A strongly deformed shock wave triggers mixing 
instabilities (Rayleigh-Taylor as well as Richtmyer- 
Meshkov) at the interfaces of the different composition 
shells of the exploding star after the passage of the out- 
going shock wave. This was shown to lead to large-scale 
mixing of the chemical elements between the deep in- 
terior and the outer stellar layers during the explosion, 
explaining self-consistently a variety of properties ob- 
served in well-monitored supernovae like the famous Su- 
pernova 1987 A, for example the high nickel velocities, 
the inhomogeneous and clumpy distribution of the met- 
als, and the spreading of hydrogen over a wide range 
of velocities that had to be invoked for explaining the 



smoothness and broadness of the lightcurve peak ||35|] . 
Big explosion asymmetries can therefore not be inter- 
preted as a signature of MHD-driven ("jet-driven") ex- 
plosions. 

Scheck et al. JsJ], performing a large set of 2D sim- 
ulations for neutrino-driven explosions with an approx- 
imative description of the neutrino transport and using 
the neutrino luminosities from the contracting and cool- 
ing neutron star as a parametric boundary condition, 
showed that such large explosion asymmetries can leave 
the compact remnant with recoil velocities sufficiently 
large to explain the measured eigenvelocities of young 
pulsars. In cases where a dipolar asymmetry became 
dominant and the explosion developed more strength 
in one hemisphere than in the other, typical kick ve- 
locities around 500km/s were found, with peak values 
even above lOOOkm/s, whereas in cases where the higher 
modes were stronger than the / = 1 asymmetry, the veloc- 
ities stayed fairly modest, usually below about 200km/s 
(see Fig. 20 in |34]). 

The pulsar recoil is caused by the asymmetry which 
the SASI distorted explosion develops on the long run, 
i.e. over a timescale of many seconds. During the on- 



set of the explosion, even until the shock reaches a ra- 
dius of some 1000 km, the pulsar kick usually does not 
grow to large values, which indicates that the ejecta have 
not obtained a high momentum asymmetry until then. 
In Fig. [6l which displays one of the cases computed in 
Ref. [34], the left panels show that the neutron star is es- 
sentially not accelerated until 250 ms after bounce when 
the maximum shock radius in this simulation is beyond 
1500km, and even at 400 ms post bounce the neutron 
star has attained a velocity of less than lOOkm/s. Only 
later the acceleration grows and leads to a recoil veloc- 
ity that asymptotes much after one second post bounce. 
The reason for such a long-lasting neutron star propul- 
sion can neither be anisotropic accretion nor anisotropic 
mass ejection in the neutrino-driven wind. The former 
ceases at about 0.5 s p. b., while the neutrino wind is es- 
sentially spherically symmetric, corresponding to a neu- 
trino emission that is nearly isotropic and thus also pro- 
duces only a very small effect on the neutron star kick 
velocity (see left lower panel of Fig.[6l). 

A careful analysis of all effects that can transfer mo- 
mentum between the surrounding gas and the compact 
remnant, i.e., gas outflow and accetion, anisotropic pres- 
sure, neutrino emission, and gravitational forces, shows 
that mostly the last mediate the speed-up of the neu- 
tron star at ^ > 0.4 s, whereas before that time all indi- 
vidual forces are large but nearly balanced (see lower 
right panel of Fig. [6] and for details, Ref. |34]). The ex- 
plosion asymmetries, which are the cause of the gravi- 
tational momentum transfer to the neutron star, develop 
only gradually when the outward moving shock encom- 
passes more and more matter from the progenitor star. 
Since the shock was launched highly aspherically as a re- 
sult of the SASI motions, the gas swept up by the shock 
may not experience the same acceleration everywhere. If 
the shock is weaker in a certain direction or is oblique 
to the radius vector, the swept up gas is less strongly ac- 
celerated. This slower gas begins to lag behind the ejecta 
expanding in other directions. It is funneled into dense, 
downward-reaching, low-entropy filaments with the typ- 
ical mushroom-like Rayleigh-Taylor caps (these are vis- 
ible in the right upper panel of Fig. [S]). The gravitational 
attraction of such massive gas pockets, which are closer 
to the neutron star than the ejecta in other directions, ex- 
certs an anisotropic pull on the compact remnant. Re- 
versely, this pull decelerates the still expanding gas, thus 
transferring some of the ejecta momentum to the neutron 
star. In the extreme case, the gas may be gravitationally 
captured and may fall back to be accreted by the neutron 
star. The acceleration ceases and the neutron star speed 
asymptotes, when the inhomogeneous ejecta reach larger 
and larger radii and the anisotropic gravitational forces 
diminish. As a consequence of this acceleration, the neu- 
tron star motion is predicted to be opposite to the main 
momentum of the matter ejected in the supernova blast. 
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FIGURE 7. Snapshots of density, absolute value of the ra- 
dial velocity, and entropy per nucleon (from top to bottom) at 
~0.56 s after supernova shock formation in one of the explosion 
simulations of a \5Mq studied in Ref. |34]. The anisotropic 
density distribution of the ejecta is visible in the upper plot, 
with mushroom-like Rayleigh-Taylor structures reaching down 
towards the neutron star at the grid center. The middle panel 
shows the absolute values of the radial velocity, which reveals 
strong sonic activity due to the impact of the downflows on 
the neutron star surface, and the lower plot demonstrates that 
the outgoing sound waves do not dissipate their energy effi- 
ciently in the rapidly expanding neutrino-heated gas so that the 
entropies of these ejecta are hardly affected. 



Alternative explosion mechanisms 



Recently, Burrows et al. 11211 12211 came up with the sug- 
gestion that supernovae might be energized by a strong 
flux of acoustic power originating from the neutron star. 
In their 2D simulations they found that the compact rem- 
nant is instigated to large-amplitude bipolar oscillations, 
/ = 1 core gravity modes, by the anisotropic accretion 
of gas. The pulsating compact remnant sends pressure 
waves into its environment, which carry a sizable energy 
flux and can even steepen into shocks, thus dissipating 
their energy in the surrounding medium and raising the 
entropy there. The ringing neutron star acts as a trans- 





FIGURE 8. Close-up of the neutron star vicinity for two 
post-bounce times (0.375 s and 0.895 s p.b.) in a simulation 
like the one shown in Fig. [7] The color coding represents the 
absolute value of the radial velocity as in the middle panel of 
the latter figure. The sonic activity by waves emerging from the 
impact of the accretion flow on the neutron star surface can be 
particularly well seen in this quantity. 



ducer that converts some part of the gravitational bind- 
ing energy released by the accreted gas into sonic power. 
The radiated sound was found to be the crucial supply 
of the developing explosion with energy and momentum 
and thus triggers acoustically driven explosions. 

While this appears as an interesting alternative to initi- 
ate the explosion if neutrino heating fails, the question is 
how the acoustic energy input compares to neutrino heat- 
ing. Burrows et al. |21] reported that in their simulations 
the acoustic energy flux dominates the neutrino energy 
deposition later than several 100 ms after bounce. Al- 
though we do not observe the large / = 1 core g-modes in 
the neutron star seen by them at late post-bounce times, 
and on the basis of our present simulations we can nei- 
ther judge nor exclude this possibility, also in our models 
— in the full-scale supernova calculations for ^S-ISM© 
stars as well as in the parametric explosions — we can 
clearly identify the presence of strong turbulence in the 
neutron star surface layer. The gas there is stirred by the 
violent impact of accretion downflows, thus creating vig- 
orous sonic activity around the neutron star (Figs. [7] and 
[8]). Very approximately, we can estimate the outgoing 
flux of sonic energy (making the assumptions of spher- 
ical symmetry and negligible energy dissipation in the 
flow) as f36l1 

^sound ^ 4nr^pv^Cs - 4nr^ (^^^ pel 

^ 5xl0^o^(^y, (3) 

where v is the fluid velocity in the sound waves, is the 
sound speed, and p'/p denotes the amplitude of the rip- 
ples on the background density p caused by the sound 
waves. The numerical value in the last expression turns 
out to be fairly independent of radius and time at some 



distance from the neutron star in Figs. [7] and [H For p'/ p 
of order unity our estimate agrees with values quoted in 
Burrows et al. f25], which might account for a sizable 
contribution to the energy of the developing explosion. 
In our simulations, we indeed see that the perturbations 
can reach such amplitudes, in particular when the down- 
flow becomes unstable and the perturbations it gener- 
ates in and around the neutron star surface layer grow. 
Outside of these transient periods, however, p'/p is usu- 
ally significantly smaller. Moreover, the net neutrino en- 
ergy deposition rates, integrated over the volume of the 
gain layer, in all of our simulations are typically several 
Iq51 gj,g/g ^ggg Eq.O, so that the energy input to the ex- 
plosions in our models is clearly dominated by neutrinos. 

Magnetohydrodynamically driven explosions have at- 
tracted renewed and increasing interest over the past 
years for a number of reasons One reason is the discovery 
of magnetars, which is taken as an indication that at least 
in some fraction of all cases very strong surface magnetic 
fields can be generated in neutron stars, possibly already 
during the evolution phases in which a supernova explo- 
sion is triggered. Another reason is the established link of 
at least some long-duration gamma-ray bursts with very 
energetic and highly aymmetric hypernova explosions of 
presumably rapidly rotating massive stars. And a third 
reason is the interpretation of asymmetries observed in 
supernova explosions and supernova remnants as conse- 
quences and relics of jet-like eruptions or as hints to the 
ejection of highly collimated material. 

Independent of how magnetic fields transfer energy 
to the explosion, e.g., by hoop stresses, magnetic pres- 
sure, or viscous dissipation of energy, the initial magnetic 
fields in the stellar core must have been amplified dur- 
ing the collapse either by magnetic field wrapping or the 
magnetorotational instability. Thus they tap the free en- 
ergy that can be stored in the highly differential rotation 
of a spinning stellar core that collapses at angular mo- 
mentum conserving conditions. This free energy, how- 
ever, is typically a rather small fraction (of order 10%) 
of the total rotation energy, i.e., 

£frf < Erot « 2 X 10^2 erg M„,,i.57?2^^ (^l^J , (4) 

where 7?ns,6 denotes the neutron star radius in 10^ cm 
and Pns is the neutron star spin period. A millisecond 
neutron star typically requires the pre-collapse stellar 
core to rotate with a period of Pini ~ ^ns(^ini/^ns)^ ~ 
10 s (Pns/ 1 ms). This is significantly faster than predicted 
by current evolution models for rotating stars, in which 
the enhanced angular momentum loss due to angular mo- 
mentum transport by magnetic field effects is taken into 
account |42]. At the onset of core collapse, the stellar 
cores are estimated to have spin periods of flni > 100 s 
(i.e., ^^ini < 0.05rad/s). This will lead to neutron stars 



with Pns > 10 ms, which is much too slow for rotation to 
be an energy reservoir of MHD-driven supernovae. Such 
a conclusion was also reached on the basis of detailed 
simulations in Refs. |38, 23]. Nevertheless, the MHD 
mechansim may still need to be invoked for explaining 
the enormous energy output of long-duration gamma- 
ray bursts and associated hypernova explosions, which 
would imply that these events are linked to rare cases 
where massive stars have achieved to retain a large an- 
gular momentum at the time when they end their lives as 
collapsars ll43ll . 

SUMMARY AND OUTLOOK 

We have reviewed a variety of recent findings that shed 
new light on the processes that cause the explosions of 
massive stars and play a role during the birth of neutron 
stars. All recent 2D simulations that were performed with 
a full 180 degree grid agree that the standing accretion 
shock becomes unstable to low-mode, nonradial defor- 
mation. This SASI phenomenon plays a very important 
role in the supernova core. It was not only found to in- 
duce a large asymmetry of the developing blast but also 
to facilitate neutrino-driven explosions by stretching the 
time accreted matter can stay in the gain layer and can be 
exposed to neutrino heating. Moreover, the SASI was ob- 
served to excite large-amplitude core g-modes in the neu- 
tron star, whose sonic damping could contribute to or be 
essential for powering the explosion. The asymmetries 
imprinted on the explosion by the SASI may lateron lead 
to neutron star kicks and might explain the observed ve- 
locities of young pulsars. In three dimensions the m 7^ 
modes of the SASI can also have an influence on the spin 
of the forming neutron star [44] . 

Despite the general agreement about the importance 
of the SASI, the physical mechanism behind this phe- 
nomenon is still a matter of vivid debate (see 1I27I1 and 
references therein) and the present calculations differ in 
many conclusions. Partly this is may be so because of 
the significantly different numerical approaches taken, 
e.g., regarding the hydrodynamics and neutrino trans- 
port, the computational grid, the description of gravity, 
and the employed equation of state and neutrino interac- 
tions. Some of the current discrepancies and controver- 
sies, however, will find an explanation when more sim- 
ulations become available and comparisons are made. 
Ultimately, however, most of the processes and conse- 
quences mentioned above will have to be addressed by 
3D models, which are currently still out of reach because 
of the enormous computational demands of the energy- 
dependent neutrino transport. 
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